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Abstract. We establish a general relation between dispersion forces. First, based on QED in 
causal media, leading-order perturbation theory is used to express both the single-atom Casimir- 
Polder and the two-atom van der Waals potentials in terms of the atomic polarizabilities and 
the Green tensor for the body-assisted electromagnetic field. Endowed with this geometry- 
independent framework, we then employ the Born expansion of the Green tensor together with 
the Clausius-Mosotti relation to prove that the macroscopic Casimir-Polder potential of an atom 
in the presence of dielectric bodies is due to an infinite sum of its microscopic many-atom van 
der Waals interactions with the atoms comprising the bodies. This theorem holds for inhomoge- 
neous, dispersing, and absorbing bodies of arbitrary shapes and arbitrary atomic composition on 
an arbitrary background of additional magnetodielectric bodies. 



1. Introduction 

The relation between the Casimir-Polder (CP) interaction of a single atom with 
a body [H [2] and its van der Waals (vdW) interaction with the atoms comprising 
the body [21 E] was first discussed in detail by Renne with special emphasis on 
an atom interacting with a dielectric half space. Approaching the problem from 
the microscopic side and modelling all atoms by identical harmonic oscillators, he 
showed that the sum of all many-atom vdW interactions between the single atom 
and the body atoms corresponds to the result for the CP potential derived earlier 
on the basis of a macroscopic approach Milonni and Lerner generalised this 
result to bodies of arbitrary shapes [6]. Using the Ewald-Oseen extinction theorem 
which follows from the Clausius-Mosotti relation, they demonstrated that the CP 
potential of an atom in the presence of a nonabsorbing homogeneous dielectric 
body can be obtained by summing over an infinite series of many-atom vdW 
potentials. 

In this paper, we approach the problem from the macroscopic side and under 
more general conditions. To that end, based on macroscopic QED in linear, causal 
media (Sec. [2]]), we first consider the CP interaction of an atom with a macroscopic 
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body as well as the vdW interaction of two atoms (Sec. [3?]), and then make use of 
the Born expansion and the Clausius-Mosotti relation to establish the microscopic 
origin of the CP potential (Sec. followed by a summary (Sec. [5?]l . 



2. Atom-field interactions 

According to the multipolar scheme, the Hamiltonian for a system of neutral 
atoms and / or molecules - briefly referred to as atoms in the following - interacting 
with the electromagnetic field in the presence of dispersing and absorbing bodies 
is given by [3 El El [TO] 

H = Y,Ha + Hf + Y.^af, (1) 

A A 

where the Hamiltonian 

HA = Y,EA\nA){nA\ (2) 

n 

governs the (unperturbed) internal dynamics of atom A, 

/f'OO 
dV / du;hu;hKr,Lo)-h{r,i^) (3) 

A=e,m ^ 

is the Hamiltonian of the system composed of the electromagnetic field and the 
magnetodielectric medium including dissipative interactions, with ix{r,uj) and 
f|(r, cj) being the dynamical variables of the system, which satisfy Bosonic com- 
mutation relations, and the atom- field interaction in electric dipole approximation 
is given by 

Haf = -dA-E(r^) (4) 
(d^: atomic electric dipole moment, r^: centre-of-mass position). In Eq. dH), 

/f'OO 
dVy da;G'A(r,r',a;)-fA(r',a;)+H.c. (5) 

A=e,m ^ 

is the medium-assisted electric field expressed in terms of the dynamical variables, 
with the quantities 



Ge{r,r',uj) = i'^G{r,r',uj) - Ime(r', cj) , (6) 
^ / / N ■ ^/ / ^ ^ / ^ Im^(r',u;) 
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being given in terms of the classical Green tensor G{r,r' ,uj), 

r.2 



Vx^ ^(r,a;)Vx — ^e(r,a;) 



G(r,r',a;)=d(r-r'). (8) 



The (macroscopic) permittivity e{r,uj) and permittivity /i(r,a;) satisfy the Kra- 
mers-Kronig relations and the conditions Ime(r, cj) > and Im /i(r, uj)>0 imposed 
for absorbing media. Note that the Green tensor has the useful properties [8j 

G*(r,r',w) = G(r,r',-a;*), G(r, r', w) = ^^(r', r, a;), (9) 
JdhGx{r,s,uj)-Gf{r',s,uj) = ^ujHmG{r,r',Lo). (10) 

A=e,m 

3. Dispersion forces 

Dispersion forces can be derived from the associated potentials U, which are 
commonly identified as the position-dependent part of the ground-state energy 
shift AE induced by the atom-field coupling. 

3.1. The Casimir-Polder force 

For a single atom A the ground state of the system is given by [0) = IOa)|{0}) 
[\0a), atomic ground state; fxi{r,uj)\{0}) = 0], and the leading, second-order energy 
shift reads 



A2E = 



{0\Haf\I){I\Haf\0) 



J Eq — El 



:ii) 



Note that the primed sum includes (principal- value) integrals over the continuous 
degrees of freedom. Inspection of Eq. ([4]) reveals that only intermediate states 
of the type |/) = |A;^)/Jj(r,Lj)|{0}) contribute. Substituting the respective matrix 
elements 

(OA!({0}|-d^-E(rA)/L(r,^)|{0})|A:^) = -[dOf.GA(rA,r,^)]^ (12) 

[d^^ = (0^|dyi|A;A)] and energy denominators Eq — Ej = —h{oj\ + oj) = {E\ — 
E^^/fi\ into Eq. pT]) . using Eq. (fim) . and separating the Green tensor into its bulk 
and scattering parts according to 

G{rA,VA,i^) = G^^\vA,VA,u:) + G^^\vA,VA,u:), (13) 

one arrives, after some calculation, at [A2E^Ua[ya)\ [2] 

UA{rA) = ^l duv?aA{iu)TTG^^\TA,TCA,iu), (14) 
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where, for simplicity, an isotropic atomic ground-state polarizability is considered, 

9 . , M I j0n|2 

= li-o 4 E , 1' . • (15) 

The potential ([H]) implies the CP force (V^ = Vr^) 

= -VAUAirA) (16) 

on a ground-state atom A due to the presence of an arbitrary arrangement of 
dispersing and absorbing bodies [which are accounted for by G(^)(r^, r^, m)]. 

To illustrate the use of Eq. (|14p . consider an atom placed at distance rA from 
the centre of a small, homogeneous, magnetodielectric sphere of radius R <tirA- 
Substituting the respective Green tensor [llj into Eq. and retaining only the 
leading-order term in R/va (cf. Ref. [E]), one obtains 



UA{rA) = - 3 2 6 / duaAiiu)[gee{urA/c)aQ{iu) 

- {uTAf /c^gem{urA/c)aQ{iu)\, (17) 
gee{x) = 2e"^''(3 + Qx + + 2x^ + x^), (18) 
gem{x)=2e-^%l + 2x + x^), (19) 

with the electric and magnetic polarizabilities of the sphere being given by [13] 

Q;?,^(a;) = 47reoi? , am('^) = ^-^^ • (20) 

®^ ^ e(a;) + 2' ®^ ^ /io /x(cj) + 2 ^ ^ 

In the nonretarded limit, where the atom-sphere separation is much larger than 
the characteristic transition wavelengths of the atom and the sphere medium, one 
may approximate geeiurA/c) ~ fifee(O), 5em(wA/c) ~ fi'em(O), so Eq. ([H]) reduces 
to 

UA{rA) = - -1^ 3 2 6 / duaA{iu)aQ{iu) (21) 

and 

UA{rA) = TTT^^ Tdn (-)'a^(m)a-(m) (22) 

for electric and magnetic spheres, respectively, while in the opposite retarded limit 
the approximations aA{iu) ~a^(0), aQ{iu) ~a0(O), aQ{iu) ~a0(O) lead to 

, , /ic«A(0)[23a|(0)-7ag(0)/c^] 
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3.2. The two- atom van der Waals force 

To calculate the vdW interaction of two atoms A and B in the presence of 
dispersing and absorbing bodies, we start from the ground state |0) = |OaO_b)|{0}) 
and consider those contributions to the energy shift that depend on the positions 
of both atoms. The leading contributions are hence contained in the fourth-order 
perturbative shift 

^E= (0|-^int|/)ai^int|//)(//[i^int|m)(///|i^int|0) 

' {Eo-Ei){Eo-Ejj){Eo-Eni) ' ^ ' 

where ifjnt = Haf + Hbf- A typical set of possible intermediate states is given by 

= \kAOB)fUr,u;)m), 1/1/(1)) = \OAiB)fUA^"',u;"')m), 

|//(i)) = |OAOB)^/l,,,(r',..')/l"."(r",t^")l{0}). (25) 
Upon using Eq. (fT2]l as well as 

- {kAOB I ({0} lA^ (r, u;)d^ •E(r^) [OaOb)/1,,, (r', u;')/l„,„ (r", w") |{0}) 
= -[[df-Gy{TA,r',oj')].,6xx"Sii''S{r-r")6{oj-uj") 

+ [df-Gx"{rA,r",oj")].„6xySu'6{T-r')6{oj-oj')} (26) 



and recalling Eq. (jlOp. substitution of the intermediate states ()25p into Eq. ([2 
leads, after some calculation, to 



2 /"oo c 

A^I)e = -^Y11 du;u;2d°fjmG(r^,r^,^).d°J 



A I '2 



d^^-lmG{rA,rB,uj')-d% ^ df ■lmGirA,rB,^')-d% 



(27) 



where we have assumed real dipole-matrix elements. Under this assumption, the 
various two- atom contributions A^^ E to the energy shift A^E only differ by the 



denominators in the square brackets of Eq. (j27p . and hence after a lengthy calcu- 
lation one may show that [A4i? i-^ [/^^^(ryi, r^)] [TOl IT4| 

UABirA,rB) = — ^/ duu'^aA{iu)aB{iu)Tt[G{rA,rB,iu}G{rB,rA,iu)]. (28) 
Jo 

Prom the two-atom potential (|28p one can calculate the vdW force on atom A{B) 
due to atom B{A) in the presence of arbitrary dispersing and absorbing magne- 
todielectric bodies according to 

Fa{b) = -y4(B)?7AB(rA,rB). (29) 
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As a simple example, consider two atoms embedded in bulk magnetodielectric 
material. Substitution of the respective Green tensor [7] into Eq. (|28p leads to 

= ~ 2.2.^elr% Jo <i^"^M«i.M , (30) 

rAB ^ jr^ — r^l, n{uj) = \J e{Lo)n{uj), recall Eq. (fTSl) . which reduces to 

aA{m)aB{iu) 

jj, N _ 23ncaA(0)aB(0) 

^™^-"64vr3e2e2(o)n(0)r7^ ^^^^ 

for nonretarded and retarded interatomic separations, respectively. Equations 
pop - (|32p show that the presence of a medium leads to a reduction of the poten- 
tial w.r.t. its well-known free-space value [2], while comparison of Eqs. (I17p and 
(jSOp reveals that in free space the dispersion interaction of an atom with a small 
dielectric sphere has the same form as that of two atoms. 



4. Relation between dispersion forces 

We now turn to the question how the CP interaction of a single atom with 
dielectric bodies can be related to its many-atom vdW interactions with the atoms 
comprising the bodies. For simplicity, we will speak of a single dielectric body in 
the following. We assume the dielectric body to be given by x(i','^)i and we allow 
for the presence of an arbitrary magnetodielectric background of additional bodies 
characterised by e(r,u;) and ^(r,u;), such that 

e(r,u;) = e(r,u;) +x(r,u;). (33) 

The Green tensor corresponding to this scenario can formally be written as a Born 
series [15j 




X G(r,si,u;)-G(si,S2,u;)- • •G(sfe,r ,cj), (34) 

where G is the Green tensor corresponding to the magnetodielectric background. 
Substituting Eq. into Eq. (fH|) . the CP potential can be written as 



oo 

Ua{va) = Ua{ya) + AkUA{rA), 

k=l 



(35) 
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where 



^vr ./n 



(36) 



is the CP potential due to the magnetodielectric background, which is not of 
further interest here, and 



AkUAirA) 



.i=l 



lUj 



xTr[G(r^,si,in)-G(si,S2,m)- • •G(sfc, ta, iu)] (37) 

is the contribution to the potential that is of A;th order in x- 

Assuming the dielectric body described by x to be comprised of atoms of po- 
larizabilities asiuj) and number densities ns(r), the gap between the macroscopic 
and microscopic descriptions can be bridged by means of the Clausius-Mosotti law 

J2B"'B{r)aB{(^)/eo 



X(r,w) 



(38) 



1 - 'Ec'^c{r)ac{uj)/{3eo) ' 

Note that since x is the Fourier transform of a linear response function, it must 
satisfy the condition x(r, 0)> x(r, iu)>0 for n>0, which implies that the inequality 



^^nB(r)aB(0) < 1 



(39) 



must hold. Substituting Eq. ([38|) into Eq. ([37|) and splitting off the singular part 
of the Green tensor according to 



one obtains 

where 

^iUAiTA) 

with 



1 / c \ 2 

G(r, r', in) = - 1 — I S(r — r') + H(r, r' ,iu) 
3 Vn/ 



AkUA{rA) = Y.^i^A{rA), 
1=1 



(40) 
(41) 



E 



m>o,-,rii>o 

rii+...+rji=k-l 



du 



X VABi...Bi{rA,si,...,si,iu) (42) 



VAi...A,{ri,...,rj)= / dnVAi...A,(ri,...,rj, 



27r 



du u'^^ 



aAiiiu) ■ ■ ■ aAj{iu)Tr[H{ri,r2,iu) ■ ■ ■ H{rj,ri,iu)] (43) 
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denotes the sum of all /-atom terms that are of order k in X; and each power of 
the factor 

,(r.i„) = - ^j.""'f^"''f;^;°> (44) 

1 - Ec^c(r)acM/(3eo) 



is due to the integration of one term containing 6{si — Sj_|_i). Summing Eq. (jU 
over k, one may rearrange the double sum as follows: 



oo k 



Y^^kUAirA) = Y,T.^'kUA{rA) = J2^'UA{rA), (45) 

k=l k=l 1=1 1=1 



A'C/(r^) 



n / ^^^^ Yi "^«(^^) 



.1=1 B, 



VABi...Bi{rA,si,...,si), (46) 



where we have performed the geometric sums X^^o means of Eq 

Note that the convergence of these sums requires \q{r,iu) \ < 1, which by means of 
Eq. ([32D is equivalent to 

^Y,nB{r)aB{0) <l. (47) 

Finally, we symmetrize the many-atom terms by introducing the symmetriza- 
tion operator 

5/(ri,...,rj) = — -i— /(r7r(i),---,Mi))' (^8) 

where -P(j) denotes the permutation group of the numbers From the 

cyclic property of the trace together with Eq. Q it follows that 

Tr [H{yi , ra , w) • • • H{Tj , ri , uj)] = Tr [H'(r^(i) , r^(2) , w) • • • H{r^^j), r^(i) , iv)] , 

(49) 

if vr is either a cyclic permutation or the reverse of a cyclic permutation. Thus the 
sum on the r.h.s. of Eq. (148^ contains classes of (2 — S2j)j terms that give the same 
result. By forming a set P{j)'^ containing exactly one representative of each 
class, we can remove this redundancy, leading to 

cSTr[i3"(ri,r2,tj) • • • H{vj,ri,Lo)\ 

= ^ Tr [i? (r^(i) , r^(2) , w) • • • -H'(r^(j) , r^(i) , w)] . (50) 

7reP(j) 

Introducing the factor 1//! in Eq. (j46p and summing over all l\ possible ways of 
renaming the variables Sj and Si, the representative of each class in Eq. ()50p is 
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generated exactly twice (only once for j = 1 + 1 = 2), so that the CP potential of 
atom A due to the dielectric body x can be written as 



UA{rA) 



oo 

1=1 



n 

i=l 



d^s, 



Bi 



UABi...Bi{rA,Sl, . . . ,Si), 



(51) 



where 



UA,...A,{ri,...,rj) = , 
{l + d2j)7r Jo 



X cSTr [H(ri , rs , iu) • • • U{rj , r i , m)] (52) 



is nothing but the j-atom vdW potential on an arbitrary magnetodielectric back- 
ground e{v,uj), fj,{r,uj). 

We have hence proved that the CP interaction of an atom with a macroscopic 
dielectric body which is described within the framework of macroscopic QED is 
the result of all possible microscopic many-atom vdW interactions between the 
atom under consideration and the atoms forming the body, provided that the sus- 
ceptibility is of Clausius-Mosotti type ([38]) and the convergence condition (|17|) 
holds - generalising the result in Ref. [6j. Conversely, our proof shows that when 
Eqs. ([5T]) and ([52]) hold and the convergence condition ([^7|) is satisfied, then the 
electric susceptibility must have the Clausius-Mosotti form. In addition, the proof 
has delivered the general many-atom vdW potentials (j52p on an arbitrary magne- 
todielectric background. For j = 2, Eq. (j52p agrees with the two-atom potential 
()28p derived in Sec. 13.2.1 while for higher j, it presents a generalisation of the 
free-space vdW potentials derived in Ref. [16]. 

The applicability of the microscopic expansion (I5ip depends crucially on the 
validity of the convergence condition (j47p . Recalling Eq. (|15p and estimating 



E 

B 



nB{r 



B 



3heo 



1 



to 



B 



/ 



(53) 



[qe, electron charge; rrig, electron mass; = h/ {00171^0); oq = q'^/{A7reQhc); -Er = 
/ {2mea\); Vs, volume accessible per atom within the body; Va = 47rag/3, / > 1, 
species-dependent factor], Eq. (f^7|) can be reformulated as ^ 8/Va/3, stating 
simply that the atoms must be well-separated within the body. If this is not the 
case, the microscopic expansion ()5ip does not converge, while the more general 
macroscopic expression (|14p for the CP potential remains valid. 



5. Summary 



We have demonstrated that on the basis of macroscopic QED in linear, causal 
media, leading-order perturbation theory can be employed to derive general ex- 
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pressions for both the single-atom CP potential and the two-atom vdW potential 
in the presence of an arbitrary arrangement of magnctodielcctric bodies. 

Moreover, starting from this very general, geometry-independent basis, we have 
used the Born expansion of the Green tensor together with the Clausius-Mosotti 
law to prove that the CP interaction of a single atom with inhomogeneous, dis- 
persing and absorbing dielectric bodies in the presence of an arbitrary magnetodi- 
electric background can be written as a sum of many-atom vdW potentials. The 
proof demonstrates the equivalence of the microscopic and macroscopic descrip- 
tions provided that the microscopic picture is applicable, while at the same time 
delivering explicit expressions for the general many-atom vdW potentials in the 
presence of magnetodielectric media. 
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